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Abstract. We prove optimal regularity for the double obstacle 
problem when obstacles are given by solutions to Hamilton-Jacobi 
equations that are not C^. 

When the Hamilton-Jacobi equation is not then the standard 
Bernstein technique fails and we loose the usual semi-concavity 
estimates. Using a non-homogeneous scaling (different speed in 
different directions) we develop a new pointwise regularity theory 
for Hamilton-Jacobi equations at points where the solution touches 
the obstacle. 

A consequence of our result is that -solutions to the Hamilton- 
Jacobi equation 

±|V;z - fl(x)|2 = ±1 in Bi, h = f on dBi, 

are in fact C^'"^^ provided that a e C". This result is optimal and 
to the authors' best knowledge new. 
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1. Introduction. 

1.1. Background. The classical torsion problem modelling the elastic- 
plastic torsion of a bar can be formulated as follows: 



Minimize dx 



in the convex set K = {u & Wj'^(Q); \Vu\^ < 1} (see [?], [?], [?], 
[?]). Brezis and Sibony [?] showed that this problem is equivalent to 
minimizing the above energy in the set K = {u E Wq'^(Q); u < d{x)] 
where d{x) is the distance function to the boundary. 

More generally, one can show that minimizing the Dirichlet energy 
in L = {w G W^'^(Q); |Vup < 1 and w = f ondD] is equivalent to 
minimizing the Dirichlet energy in the set L = {m G W^''^(Q); h~ < u < 
h'^} where h- solves the Hamilton-Jacobi equations 

±|V/z±P = ±l inQ, 
h^=f on dn, 

provided that L and L are non empty, which incidentally are equiva- 
lent conditions (see [?] for a proof of this equivalence). 

The regularity theory for a minimizer of the Dirichlet energy in 
L is quite straightforward. Indeed, one may approximate h'^ by the 
solution /i+ to 

-eA/i+ + IV/i+p = 1 inQ, 
h+ = f on dn, 

whence the Bernstein technique (see [?] or Lemma \&2\ below) gives 
an e-independent estimate on the second derivatives from above. In 
particular, letting e — > we may deduce that h'^ is semi-concave, 
or equivalently that the distributional second derivatives of h'^ are 
bounded from above. Similarly, the second order distributional 
derivatives of h~ are bounded from below. From here, standard 
regularity theory for the obstacle problem (such as developed in [?], 
modified slightly in the present paper to suit our purposes) implies 
thatw G 

loc 

The important step in the above proof is to deduce one-sided es- 
timates on the distributional second derivatives of h-. More gener- 
ally, if M is a minimizer to the Dirichlet energy in L with h'^ being a 
solution to a Hamilton-Jacobi equation F(x, h"^, V/i"^) = 1 {h~ solves 
-F{x, h~, V/i") = -1) for an F G satisfying some structural assump- 
tions, then the same technique yields u G C^'^. The important step 
again is the Bernstein technique where we need to differentiate F 
twice. See for instance [?] or [?] for variational problems of this 
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type or [?] with errata [?] for a fully nonlinear gradient constrained 
problem. 

WhenF ^ C^, the above outline of a regularity proof fails. The exist- 
ing regularity theory for Hamilton-Jacobi equations is not sufficiently 
strong for that purpose. There are cases where certain one-sided esti- 
mates have been deduced, see [?] and [?]. For example, in [?] and [?] 
it is assumed that f G C°'^ (besides standard structural assumptions) 
and it is shown that h'^{x + x^) < + p ■ x + |x|a(|x|) for any p in 

the super-differential of h'^ at x° where o is some one-sided modulus 
of continuity over which we have no control. 

1.2. Main Result and Ideas. The objective of this paper is to intro- 
duce new techniques handling regularity questions for Hamilton- 
Jacobi equations below the -threshold. It is noteworthy that the 
minimization problem in class L with Hamilton-Jacobi equations be- 
low C^-threshold has applications in micro-magnetics. In view of the 
equivalence of the minimization problem in class L and L for the spe- 
cial case of gradient constraint |Vm| < 1 (see Background), we study 
in the present paper the problem in L with variable coefficients that 
are not Lipschitz. 

Our main results are the following theorem and its corollary. 

Theorem 1.1. [Main Theorem] Let ubea minimizer of the Dirichlet energy 

f \Vuf 

in the set {u e W^'^; h~ < u <h'^, u = f e C"{dBi) on dB^}, where h- are 
viscosity solutions to 

, . ±|V/z± - fl(x)|2 = ±1 mBi, 

^ ' =f on dBi, 

with [^z]q,(b7) = a < C for some C < +oo. Assume furthermore that 
m(0) = h+{0). Then 

osc^eB,{o) (u{x) - m(0) - Vi/(0) • x) < C{a) y/Ar^+'i ifr < A^^, 
osCxeBr(o) {u{x) - m(0) - Vi/(0) • x^ < Cr^ ifr > A^. 

A surprising consequence which may be of immediate interest to 
the regularity theory for Hamilton Jacobi equations is the following 
corollary. 

Corollary 1.2. Let hbe a solution to \Vh - a\^ = 1. Then h e C^'"^^, 
provided that a e C^. 
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Proof. Since h e C^, we have, by uniqueness of solutions to HJ- 
equations, that if h- solves equations 

±|V/z± -fl(x)P = ±1 inBi, 
h- =h on dBi, 

then Ifi" -h~ - h. In particular the set 

K = {ue W^'^; h- <u<h*] = [h]. 

Therefore, is in a trivial way a minimizer of Jg |Vu|2 in K. By 

Theorem O it follows that h e C^'"'^. □ 
The function F{x,p) = \p - a{x)\^ is related to a HJ equation that 
arises in micro-magnetics, and hence our choice is not completely 
arbitrary (see [?]). Our method is quite robust, and as such it seems 
plausible to adapt it to a wide class of Hamilton-Jacobi equations. 

The main difficulty, as indicated above, is to develop a strong 
enough regularity theory for Hamilton-Jacobi equations. In this arti- 
cle, we will not treat Hamilton-Jacobi equations in their full general- 
ity. Instead we will stay within the confines of the obstacle problem. 
This has one great advantage: it is easy to see that u G C^'^ for some 
jS > (Lemma 16.51 and Proposition 12. 5|) . We will therefore get a 
one-sided estimate of h'^ from below at all points where u = h* . That 
means that the set {w = h+} does not intersect the singular set of h.-; by 
the singular set of /i"^ we mean the set where h- are not differentiable 
in the classical sense. 

The regularity theory for h"^ is deduced by inhomogeneous scaling. 
There is a slight complication to apply this method to Hamilton- 
Jacobi equations of our type. In particular, even though h'^ satisfies 
an elliptic Hamilton-Jacobi equation it scales parabolically and the 
blow-up limit will solve a parabolic equation. Let us denote 

h+{rjX',r]~'^x„) 

hM) = t4 ■ 

r. ' 

Then if hj — > ho as Tj — > 0, the function ho heuristically solves an 
equation of the form |Vhop - 2d„ho = where V = (di, ...,dn-i,0). 

Our first goal (Proposition I3.2|) is to show that the hj defined above 
is indeed bounded. It is here that we use the assumption u(0) = h'^{0), 
which gives a one-sided estimate from below on h'^. 

Once that is proved we can use the regularity theory for para- 
bolic Hamilton-Jacobi equations to deduce that h'^ satisfies some- 
what better one-sided estimates in the x' directions, for example 
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h'^{x',0) < /z+(0) + p ■ x' + Qx'l^'^^'^^ where e - e{(i,a) (see Proposition 
iH). 

The minimizer u on the other hand scales elliptically, so if 

u{rjx)l sup \u\ 

Br. 

I 

converges to Uq as ry — > 0, then Uq is a solution to an obstacle problem 
with obstacle Iiq = limy^oo h'^{rjx)/ supg \u\. This fact will be used in 

the proof of our main theorem. Theorem 11.11 in order to show that 
oscb,. \u{x) - Vm(0) • x\ < C\x\^'^f^'^^ . In particular we gain an e in the 
regularity of u. By carefully keeping track of all the constants we see 
that this can be iterated to obtain y-independent C^''' estimates of u 
for all y < a/2. It follows that u G O'"'^. 

The following example pointed out to us by Stefan Miiller shows 
that this is indeed the optimal regularity: 

Example 1.3. Ifa{x) = \x2\"ei, then u{xi,X2) = Xi + b{x2) is a solution to 
|Vw-fl(x)P = 1/or 



HX2) = 




is a solution to |Vw - a{x)\^ = 1. Here a e C" and u g C^'"/^, hut u ^ C^'^ 
for any jS > a/2. 

The plan of the paper is as follows. In Section |2] we deduce an 
abstract regularity result for solutions to the double obstacle problem, 
which we will need later. In the subsequent two sections we show 
that remains bounded under parabolic scaling and that h'^ satisfies 
better one-sided estimates in the x' directions. In the final section we 
prove our main result that the minimizer u G C^'"^^. Finally we 
have included a long appendix to remind the reader of some of 
the theory of viscosity solutions for Hamilton-Jacobi equations. In 
the Appendix we also deduce simple C^'/'-estimates for the solution, 
which will serve as our starting regularity in the strategy described 
above. 

Acknowledgment: We would like to thank Stefan Miiller for pro- 
viding us with Example II. 3[ 

1.3. Notation. We denote the Euclidean ball Br{x) = {y G R"; \y-x\ < 
r\, and we denote cOn '■= £'"{Bi); in the case that the center is not 
specified it is assumed to be the origin. When v e M" is a vector 
we will denote the first n - 1 coordinates by v' := {vi,V2, ...,Vn-i). 
Similarly, we will use V := {di,d2, ...d,j_i,0). Here dj = denotes 
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differentiation with respect to the x, variable. We will often denote 
differentiation by a subscript djU =: Uj. The unit vector in the z'-th 
coordinate direction will be denoted by e,. We will also use the 
seminorm 

\a{x) - a{y)\ 

[a\a = l«Jc«(Q) •= sup . _ . ^ • 

.v-,j/6Q \^ y\ 

Finally, different instances of the letter C may mean different con- 
stants even within one proof or one set of inequalities. 

2. Regularity for Double Obstacle Problems 

In this section we will prove regularity for the two-sided obstacle 
problem in the context needed later. 

It will be convenient to define a class of solutions to double obstacle 
problems and to fix some notation. We therefore state the following 
two definitions before we state and prove the main result of this 
section. 

Definition 2.1. For any continuous function u we define the super-differ- 
ential (sub-differential) ofu at the point x° G Domain{u) as follows: 

S*{u, x°) = {p G R"; sup [u{x + x°) - m(x°) - p • (x - < o(r)) 

B,. 

(s-(t/,x°) = {p G R"; inf (u(x + x°) - - p • (x - x°)) > -o(r))). 

Definition 2.2. We define G{R, a, h-) as the set of local minimizers u to the 
Dirichlet energy 

f Nuf 

Jbr 



in the set K = {u e Wj^^(Br);^" < u < here < R g R, o{r) is a 
"one-sided modulus of continuity" ofu and h- G C°(Br) such that h+ > h~. 

Furthermore we require that u satisfies the following estimates for each 
point x° G A{u,R) := ({w = h+} U {u = h'}) n Br; 

(i) lohen w(x°) = /z+(x°), then we have the bound from above 
sup inf (u{x) - p^ ■ {x - x^) - u{x°)) < o{r)for all r > 0. 

(ii) lohen m(x°) = h~{x^), then we have the bound from below 
inf sup (m(x) - pxo ■ (x - x°) - m(x°)) > -o{r)for all r > 0. 
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Remark 2.3. It is important that by (i) and (ii), u inherits its "one-sided 
modulus of continuity" from that ofh-. 

Proposition 2.4. Let u e G{R,o,h-) where o{r) is a one-sided modulus of 
continuity satisfying the doubling condition 

o{Sr) 

(3) — -— < f (S) for each S > 1 and all r > 0. 
o{r) 

for each S > 1. Then ifx^ e A{u,R/2) it follows that 

osc inf (u - Vyo ■ (x - x°) ) < C(n, F)o(r). 

Remark: It is of no importance that the set K has the particular 
form K = {u e W^'^; h~ < u < h'^}. We also do not need that h- 
solves any particular Hamilton-Jacobi equation in this Proposition. 
The important thing is that we have one-sided estimates at points 
where u is not harmonic. Similarly the boundary values / are of 
minor importance. 

If we consider the function F in (|3l) to be F(s) = then the Proposi- 
tion slightly generalizes known regularity results for double obstacle 
problems. 

Proof of Proposition \2A\ We argue by contradiction and assume 
that there is a sequence Uj of minimizers as in the proposition, with 
o = Oj {Oj corresponds to a fixed f ) and points x' G A{u, R) (without 
loss of generality we are going to assume that x' = 0), pj G S'^{u, x) or 
Pj G S~{u, 0) as well as ry > such that 

Uj{rjX) - rjPj ■ X 

(4) osc inf : — — = 1. 

We may also assume that for large /, is the largest such r corre- 
sponding to u', that is for p G S'^{u, 0) or p G S~{u, 0) 

Ujirx) - rp ■ X 

(5) osc " . ; < 1 

for r >rj. 

Next, we define 

Uj(rjX)-rjPj-x-Uj(0) 

where is the vector in equation dD. Thenu^ G Q{R,o{rj-)/{jaj{rj)),hj) 
where 

h'jlrjx) - rjPj ■ X - Uj{0) 
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In particular Au^ = m some set Qy = Br/,-. \ A{v>, R/rj), and for every 
point x° G A{v', Rl{2rj)) we have the following: if v^{x^) = hj{x^) then 
(cf. condition (i) in Definition 12 .21) 

(6) sup inf (v'{x) - p • (x - x°) - v'{x°)) < < — ^ 0. 

Moreover, if v^{x^) = h'j{x^) then it follows from condition (ii) in 
Definition |Z2] that 

(7) sup inf sup (v'{x) - p • (x - x°) - v'{x°)) > -"^J^ <^}. ^ q_ 

Let us write A{v>,R/rj) = A+{v',R/rj) U A-{v^,R/rj) with © satisfied 
for points in A"^ and (0) satisfied in A". 

For some subsequence — > in L|'^^, where 

for alH < +oo; here we have used @ and <^ which imply that 

o{rjr)/{jOj{rj)) ^0 as ; ^ oo, 

as well as the notation = limy^oo^*- We will actually show later 
that this convergence is locally uniform. Furthermore, by equation 
© we have that S+(u°,x°) 9^ for all x° G A+(i;°, +00) and 

v%x)< inf inf (v%x^) + p\-{x-x^)\ 

and by equation we have that S"(u°, x°) for all x° G A"(z;°, +00) 
and that 

v°{x) > sup sup (z;°(x°) + p^o • {x - 

That is, solves the double obstacle problem with a concave upper 
obstacle and a convex lower obstacle. It follows that Av^ = in R". 
By our assumption that the origin is in A(u°, +00), say G A'^(u°, +00), 
it follows that for p G S+(z7°, 0) 

sup - p • x] = 

Br 

for all r < +00. 

It follows from Liouville's Theorem that v° is a linear function, but 
i'°(0) = and by construction (cf. equation (ID) 

(8) osc inf (vHx) - p • x) = 1. 
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Equation ((8l) is a contradiction to being linear, provided that we 
can show that v' — > uniformly. 

In order to show that — > v° uniformly we notice that the obstacle 
functions 

/+:= inf inf (v'{x^) + v,o ■ {x - x^)) + -r-^ 

' xOeAJp^oeS^vi,xO)^ ' ' ]o{rj) 

converge either locally uniformly to +00 or that locally uniformly 
limy^oo = where is an affine linear function. Similarly, 

:= sup sup {v'{x'') + p^o-{x-x'')j- —— 

will converge either locally uniformly to -00 or locally uniformly 
limy^oo /" = /" where /" is an affine linear function. 

Since fr < u' < f+ we have f r < f+. 

We may distinguish three cases: 

1) v\x') = fj{x^) for some bounded sequence of points x^ and 
v^{y^) = fj~iy') for some bounded sequence of points y', 

2) v'{x') = fj{x^) for some bounded sequence of points x^ and 
> fj in Br, with Rj — > +00, or v'{y^) = fj{y') for some bounded 

sequence of points p and < in Br, with Rj — > +00, 

3) /r < < fj in Br, for a sequence Rj — > +00. 

In the first case it is easy to see that \ fj - /,~| — > locally uniformly: 
if not then /r — > a+p-x and fj — > b+p-xwithb > a. Alsoy^ — > c+p-x 

in L^^^ with b > c > a. We may assume that p = 0. Now only one of 
the inequalities c > a + {b - a) /3 or b - (b - a) /3 > c may hold; let us 
for definiteness assume that b - (b - a)/3 > c. Then since — > c in L^^^ 
there is for T := 2 supymax(|x^|, \y'\) and each e > a je < +00 such 
that 

(9) ^£,(„,-.rf<. 

for j > je- Also, as fj converges uniformly, fj < c + e for all / > je 

(provided that je has been chosen large enough). In particular {z;^ > 
c + e} n At =0 when / > je. That means that may touch only the 

upper obstacle in {v' > c + e}. Thus v' is subharmonic in {z;^ > c + e}. 
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From (ID we deduce that there exists f G (37/2, 2T) such that 



■J— f {v^-cy 



(10) — I (z;^ - cr < Ce 

for some fixed C > 1. Now we define zv^ by the Poisson integral as 

ir \ - T max(z;^, c + Ce) 

(11) w!{x) = ^ -— — dy. 

Then Aw' = in and lu' > v' on dBf. Since A.v' > in {i;^ > c + e} 
and w' > c + Ce it follows that max(u-', w') is subharmonic. But lu' > v' 
on dBf and Azv' = so max{v', lu') < w' < max(u^, w'), where the first 
inequality follows by comparison and the second is trivial. Thus 
w' > vL Moreover we know from (|TT|) and ((T0|) that w' < c + Ce in Bf. 
If e is small enough we may deduce that in Bf, 

^ , h - a ^ , h - a 
v' <c + Ce <b — + Ce <b — < /+. 

3 6 ^' 

But that contradicts the condition in 1) that v' touches the upper 
obstacle. 

In case 2) we may argue similarly as in case 1) to show that for each 
T, fj > v' > fj - e in Bj if ; is large enough, or fj < v' < f7 + e mBj 
if ; is large enough, respectively. 

In case 3) we have that for each R < +oo, v' is harmonic in Br for 
large /, which implies uniform convergence in C^^^. The Proposition 
follows. □ 

Using Lemma [631 we can derive a simple Corollary that states some 
preliminary regularity for solutions to the obstacle problem. 

Corollary 2.5. Let u be a minimizer to the two-sided obstacle -problem as 
above, i.e. u is the minimizer of 



/ 



I Br 

in the set [u g W^'^; h~ <u < h'^} and h- are solutions to ±\Vh- - = ±1 
where [a]a < A. Then ifx^ G Br/2 and w(x°) = /z+(x°) (or u(x°) = /z"(x°)) 
and is in the super(sub)-differential ofu at it follows that 

(12) osc {u{x + x°) - ■ (x - x°)) < c(n)Ai/(2+a)^i+«/(2+«) 
for r < CA'^'^R^ and that 

(13) osc (u{x + x°) - p^o ■ (x - x°)) < -r^ 



o{r) = 
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In particular, u e C^'^{Br/2). 

Proof: Using Lemma [631 we see that u satisfies the assumptions of 
Proposition 12.41 with o defined by 

for r < CAi/2R(2+a)/2^ 

o{r) = C-, 

for r > CA^/^R^^+"y\ 

It follows then by a standard method that u G C^'^: from (|T2)) 
it follows that Vm(x°) is well defined. Whitney's extension theorem 
implies that we can extend u\ , to a function u defined in Br/t 

^ isupp(Ai() ^^'^ 

where Vu is C" with modulus of continuity 

C(n)Ai/(2+«)r«/(2+«) for r < CA^'^R"^, 
^ §r for r>CAi/2i^^. 

In particular if x, y 6 supp(Au) then |Vw(x) - Vu{y)\ < o{\x - y\). 

Also notice that if dist(x°, supp(Aw)) = - y\ = k where y G 
supp(Au) then Aw = in B^ix^) and the supremum of \u{x) - u{y) - 
Vu{y) • (x - y)| in B,.(x°) can be estimated by o{2k). Standard regularity 
theory for harmonic functions implies that |Vi/(x°) - Vi/(y)| < Co^Ik). 
In particular it follows that if x G Br/2 \ supp(Ai/) and y G Br/2 n 
supp(Aw) then |Vt/(x) - Vw(y)| < Co{l\x - y|). 

If both y, z G Br/2 \ supp(Aw) and 
1 

\y - z\ < - max(dist(y, supp(Au)), dist(z, supp(AM))) = k, 

then u is harmonic in B^{y) or in B^{z). For the sake of definiteness 
we will assume that u is harmonic in B^iy). Also supg (^j|i/(x) - w(y) - 
Vu{y) • (x - y)| < o{2k). Standard regularity theory for harmonic 
functions implies that 

\Vu{y) - Vm(z)| < Co{2K)\y - z|/k < Ca{\y - z|). 
Finally, if y,z G Br/2 \ supp(Ai/) and 
1 

ly - z| > - max(dist(y, supp(AM)), dist(z, supp(AM))) = k, 

then we may combine the first two estimates and deduce that 

\Vu{y)-Vu{z)\<Co{\y-z\). 
The Corollary follows. □ 
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3. Improved Regularity in the Non-Degenerate Direction 

In this section we will discuss the behavior of a viscosity solution 
in the non-degenerate direction. To motivate this discussion and our 
terminology, let us consider a simple elliptic rescaling of | V/z"^ -ap = 1: 
if is differentiable at the origin, /z+(0) = |V/i+(0)| = 0, a{x) = e„ + 
\b{0)\ = and 

' ' supg l/i+l 
converges to Hq as fj — > 0. If we rescale we see that 

supe J/z+l 

=|V/iy|2 1 - Ihifjx) ■ Vhj 

+ Hrjxf ^ + IbniTjX) 

^ supg^ |n+| ' supg^ \h+\ 



and — at least heuristically- 



dx„ 



provided that 



\birjx)\' '-— + IKiTjX) '-— ^ 0,; ^ oo. 

' supg \n+\ ' supg |n+| 

'7 '7 

Therefore the Hamilton-Jacobi equation is degenerate in all directions 
except one. Thus we might expect that the oscillation of h'^{x' ,x„) 
along lines {(x',s); |s| < r] with fixed x' should be of lower order 
compared to oscg^ . In the next Proposition we will show that is 
indeed the case, if we have a good estimate on from below. 

The proof, and even the statement, is quite long and not very easy 
to read. Therefore we will try to explain the general idea before we 
start. 

The idea of the proof is that scales like h'^{rx', r^~f^x„)/r^'^f^, when- 
ever h+{0) = |V/z+(0)| = (notice that if u{0) = h+{0) then V/z+ is 
defined at the origin). It is therefore a natural assumption that 
h'^{rx',r^~f^Xn)/r^'^^ will be bounded if u G C^'^ and h'^ has one-sided 
C^'^'^ -estimates. We will prove this boundedness (or a slightly refined 
version of it) by blow-up and a contradiction argument. Most of the 
proof consists of technical estimates of different terms in the equation 
for the rescaled h'^. The idea is simple though. If the supremum of 
the rescaled h'^ goes to infinity then we can find another rescaling 
(called hj in the Proposition) that respects the natural scaling of the 
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Hamilton-Jacobi equation (Claim 1, below) and is bounded. This 
new rescaling of h'^ has the nice property that it is worse than a C^'^ 
scaling in the x' directions. This will imply that hj — > ho where ho 
is independent of the x' directions, at least at x„ = (Claim 2 in the 
Proposition). The remainder of the proof consists in verifying that ho 
also satisfies a good Hamilton-Jacobi equation that will allow us to 
conclude that = which in turn will imply a contradiction to the 
fact that ho{0) = and supg^ \ho\ = 1. 

Remark 3.1. In the sequel lue may assume that A > 0. The case A = 
may then he handled by approximating < Ay — > 0. 

We will denote 



Keeping this definition in mind, let us now state the main result in 
this section. 

Proposition 3.2. Let h'^ be a viscosity solution to |Vh+ -a\^ = 1 in Br such 
that h'^{0) = 0. Assume furthermore 

(al) that h'^ satisfies one-sided C^'^ estimates: 



for every x° and every p in the super-differential ofh'^ at x^, 
(a2) that 



(14) 



(15) 




(16) h*{x',Xn) < h*{Q,Xn)+p- x' + CA^Ix'r^ 




Then, for r < min(A^/2j^i-/',R), 



sup \h^{x) + g{x„)\ < CA^r^+^ 



m 



where 



(18) g(0) = g'(0) = 0, g'(0 = -en ■ b{0', t), 



and K{r) is as defined in (fill ). 
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Remark 3.3. The condition r < mm{A^^^R~^ ,R) assures that K{r) c Br 
but is otherwise not used in the proof. 

Remark 3.4. The assumption (a5) may always he satisfied scaling and 
rotating, so that it is not restrictive at all. 

Proof. We will argue by contradiction and assume that there exist 
sequences ry, Aj, V and a^ satisfying the hypothesis, such that 

suplh. +gj\ = ]A ' r 

K{rj) 

where gj satisfies (|T8|) . The proof is rather long so we will divide it 
into several claims. The first one will slightly modify K{rj) so that it 
respects the natural scaling of the Hamilton-Jacobi equation. 

Claim 1: For 

Sj{o) = sup {hj{ yjdrjx', r'-^A'^Xn) + gj{r]r^ A'f^ Xn)), 

hi i+fl 

there exists Oj such that Sj{oj) = OjA.^ r. ^ and Oj — > +oo. 

Proof of Claim 1: By assumption Sy(l) = jA. ^ /.^'^ » A' Set- 
ting Tj := 1/(aJ ^r^^) and using the one-sided estimates for h'^ (in- 
equality ((15))), we also obtain that 

Sj{Tj) < sup \hj\ + sup \gj\ 

^,i-/i^-(i-;i)/2 ^,i-fi^-(i-ff/2 

/ i i i 

< CA.^ r.^ +/.^ 2 

< c(x,r';^Af)Af^r^f < x.J^-^''^^ ' 

provided that r^ is small enough. 

By the continuity of Sy there is a 6 (l/Ty) such that Sy(c7y) = 

cr. A. r. ^ . 

! 1 1 

Therefore we only need to show that cry — > +oo. By definition we 
have 

hji Va-r,x', r]-^A-^'-'^'\) + 
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A;(r,a',r;-''/l-"-'»«x.)+ft(r;-'A-"-««x„) 

- T ^"^^ = 

] 1 

which proves Claim 1. 
Let us also remark that as 

SjiOj) = a^Afr^, 
we obtain from our one-sided estimates that 

Sj{pj) < Cmax( sup /z| + sup sup /z| + sup 

Notice that, since g'.{t) = b{0', t), 

sup \gj\<C\A. r\ '"jA.^r.', 



but the term in the parenthesis tends to zero, implying that 

ill i_«2 
sup \gj\ < A. ' r. 



B 

I I 



Similarly, 

sup |^,-| < Ajof^"^^V-^^ < (ofAfr"-^))ajAfrY^ < Sj{oj). 

Using these estimates on \gj\ we arrive at 

S,(c7,) < Cmaxfcj.^ A.^ r.^'^,A.^ r. ^ ). 
Since ay — > +00 it follows that 

cr.A. r.< CA. ^ r. . 

We may — in order to simplify notation — assume that Oj = j and 
rewrite the above relation as 

(19) Tj < CA. j . 
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As mentioned before. Claim 1 gives an estimate that respects the 
scaling of hj: if we define hj by 

h (x) ■= — - - - — - 

lA.^ r. ^ 

' 1 1 

then supg \hj\ = 1. Also 

where y = iy^fjrjX'j^. ^a/^ ^''^'^x^^. Alternatively we may write 
where 

Rewriting the Hamilton- Jacobi equation ([T]) in terms of our new func- 
tion hj and using assumption (a5), we may deduce that 

|Vh/ - Idnhj = 

(20) ^ + 2dnhj{h ■ en + g ) ^ - 2— - c5(a„/Z;)2 = 

In order to use this equation we need to control the right-hand side. 
First we need to control the one-sided oscillation in the x' directions. 
Claim 2: For every x„ and every pj in the super-differential of hj at 

{x', x„) we have for each R and any qj in the super-differential ofhj at 

(0,rJ"^A^^) that 

sup (hj - hj{0, Xn) - pj ■ x'j ^ as j ^ oo. 

In particular, if Xn = thenoscB'^{hj{x' ,0)-hj{0)-pj-x') as j ^ oo. 
Proof of Claim 2: The first statement is deduced from (a2) in the 
following way: 

h+ -h+{0,Xn) + cjj-x' pi+p 
sup (hj - hj{0,Xn) - Pj ■ x'] < sup — < 

b; b' , jA—y^+P j— 
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as oo. For the second part we use ([TT)): 

hUx',0) pi+f^ 
mfhj{x',0) > inf > -C^-r ^ as / ^ oo. 



'V/'yS ;A.^ r^^^^ r 



In order to prove the proposition we want to show that hj — > /zq = 0, 
which would clearly contradict uniform convergence and the fact that 
by definition supg^ \ho\ = limy supg^ \hj\ = 1. 

To show the convergence to we need to control the right-hand side 
inequation ((201). In particular, we are going to prove that r^,...,Tg — > 
as y — > oo. We formulate this as a new claim. 

Claim 3: We have — > locally uniformly as j — > oo, for i = 1, 5, 
where Tl has been defined in equation |20D . 

Proof of Claim 3: Using that \b{x)\ < A\x\" and that ^<^< a/2, we 
see that 

sup |r;i = sup — — - < — — = > 



as j — > oo, which implies the claim for T|. 
We also have, by assumption, that b{0, x„) • e„ + g'j{x„) = 0, and thus 

(21) sup \b{x) ■ en - g'j\ < sup osc (b ■ e„) < Ajf^^r" 
Using Lemma [6l6l we see that 



|^„/i++g'.| cA. ^ f. 
(22) sup|<9„/zy|< sup hnr^ ' ' 



(x'|{^fjrj),A^Xn|(rY'))'^B^ ' ' 

From (IZD, CZl) and ]6 < «/2 < 1/2 we see that 



sup|r^|<c^ 



as y — > oo. Next, we estimate in Bi 



(23) |T^|<— Ll_< 



2|&| ^ /■ /■ /■ 



'11 ■'11 
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as y — > oo. The term satisfies 

\T[\< " f/ < sup osc (^l^L-) < -L/ ^ 



as / — > oo. 

So far we have shown that Tj, as /—> oo. 

Next, we derive estimates for T^. First we use (|T9|) to arrive at 

' ^ i i 

Together with Lemma [6l6] this implies that 

in{x; {x'/{^fjrj),Afx„//-^)eB^l 



|V/zt| < CA.^ A. ^'r 



that is 
(24) 



|%| < C- 



Now, since /zy is a Lipschitz solution to the Hamilton-Jacobi equa- 
tion, the super-differential consists almost everywhere of a unique 
element p, and almost everywhere, Vhj = p solves equation (|20|) . That 
is, at almost every point. 



jA'-^'.^' jA'-^'P 
'11 '11 



0, 



where 



6 = |Vh/ - 



• Vhj - Idnhjiibix) ■ en + g'j) 



\h + engf en-b + g'. 

'11 '11 
This means that p„ can only take the two values (depending on 6), 

(25) pn 



' 1 I 



? + 



' 1 I 





6 




' 1 1 




6 




' ! 1 



NON-C^ HAMILTON-JACOBI OBSTACLES 



19 



By (|22l) as well as the above estimates on |V/zy| (c.f. equation (|24l)) we 
see that 



(26) 



|6|<C 



A 



a: 



l-a+a/i+/J^ 



^«-^(l+a+)3) 



+- 



■l_n 



+ 



+ 



■1_n 



Next, using (|T9|) and (|23ll we see that when / is large and fj is small. 



■l_n 



r 



Similarly, using (|19|) which implies that 

a+ji-l 

i~A. 



we may conclude that 



a: 



(28) 



^a-p(l+a)-^ 

i 



,2 l-a+a/i+/J^-2/J 



V y«/2 / 



-a/2 



< j i+P 2. ' ' 



1 



jl-a/2 ' 



Noticing that the term ; ^ in (|28)) goes to zero as 
deduce that 



A 



(29) 



<< 



A^.r"-'^ 

1 1 



1 



1 



l-a/2 



as oo. 



Moreover, as 2j6 < < 1, we have 

^l+l}-{l-fi)a^2(a-fi(l+a)) 



« 



(30) — ^ 

when Vj is small. Using (|26|) - (|30l) , we estimate 



|5|<C 



A 



-(l-p)p-2fi- 



1 



r. • A'r 



/3 a-2p 



■1-" 



oo, we 
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Using (|19|l again along with a/{2 + a) < (} < a/2 < 1/2, we see that 



a: 



-2/J2 



(31) 

and that 
(32) 



<< 



/ 1 



2/i-£V + It/J 



A:(^-«r:^^ 



<{— — ]- — - 

\ -a(fi+f-)+2a I j 



<< 



when i is large. From (|25l) , ((3T|) and (|32l) we see that we can make a 
Taylor expansion in equation (|25l) and deduce that p„ is of order 

2 



(33) 



\Pn\ 



' I 1 



2 - i 



(1-/J)B 



+ 



Now remember that from (|22|). we have the additional information 
that 

cA. ^ r.^^ 

\dnhj\ < . 

; 

As j6 < 1 /2, we conclude that for small rj, the first line in (|33)) and not 
the second must hold. Thus 

|p„| = + lower order terms. 



It follows that almost everywhere, 

T3 = jA]-^fip„f < C 



i i , /■ 



+ 



l-a 



Since 2(} < a < 1, the right-hand goes to zero as / — > oo, and the claim 
follows. 

The proof of the proposition is now easy to finish. We know from 
Claim 3 and Lemma [531 that hj — > Hq locally uniformly in R", where 
ho is a viscosity solution to 

(34) - 2d„ho = 0. 

Moreover hj — > on {x„ - 0} by Claim 2. Using uniqueness for 
parabolic Hamilton-Jacobi equations ([?]) we deduce that = for 
Xn < and /zq > in {x„ > 0}. By Claim 2 we also know that h^ 
is concave in the x' directions, which together with the bound from 
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below implies that ho is constant for each x„ > 0: if ho{x',x„) is not 
constant for some x„ = t > 0, then ho{x', t) < ho{0, t) + p ■ x' for 
some non-zero p G R""^. It follows that there is a point such that 
hoix^, t) < 0. But from (|34ll it follows that d„ho > 0, a contradiction. 
Therefore ^o(^) = hQ{x„). Using ((34|l again we see that d„ho = in 
{x„ > 0}, and thus ho = 0. This contradicts supg^ \ho\ = 1, and the 
Proposition follows. □ 



4. Improved Regularity in Degenerate Directions 

In this section we improve the regularity in the directions orthog- 
onal to the non-degenerate direction. The idea is to use the scaling 
from Proposition 13.21 to get a parabolic equation and deduce better 
regularity in the x' directions from the regularity theory for parabolic 
Hamilton- Jacobi equations. 

Proposition 4.1. Let h'^ be a viscosity solution to |V/z+ -af- = \ in Br and 
assume that h'^ satisfies the one-sided C^'^ estimate from above 

(35) sup (/z+(x) - h+(0) - p • (x - 0)) < CA^ri+^. 

B,(0) 

Furthermore let 

(36) h*{x',Xn) < h+{0,Xn)+p- x' + CA^Ixf+f* 

for each fixed x„, and assume that h'^ satisfies the one-sided C^'^ -estimate 
from below at the origin 

(37) /i+(x) > -CA^|xp+^. 
Also assume that 

(38) fz+(x',0) > /z+(0) + p-x' - CA^Ix'r^ 

In the just mentioned estimates we assume moreover that the constant 
A < C (where C is a fixed constant) controls the seminorm of a{x) by 
[fl]c« < A. Finally we assume that a{x) = e,, + b{x) where b{0) = and that 
a/{2 + a) <^< a/2. 

Then for every every 6 < min{A^^^R^ ,R) and every x„ G {Q,A~^6^~^) 
we have for some p g R" 

sup ih\x)-h-'{Q,Xn)-p-x') < Q 
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where 

CA— 6^+""^ z/A ^ < 6 < A , 

C^i+zi^i+i+p z/A < 6 < min(A ^/i ,A— R^+^'j, 

CA— 6^+^^ z/min max(^A 2/'', A w-a) 



Q = 



6/i+a/i-2-a v 



< 6 



ifR = mm{A^'^R^P,R), and 

( 1 1 , a / 6/j+a/i-2-a g_ 2+a ?:L£Zr£\ 

^^2 2(2+a-2/3)5i+2+a-2/3 § < mill ( A ^W""' ,A «-2/i , A" ^(l-fl 

2^ 2+a-2/i ep+afi-l-a 

z/A— < 6 < A , 
Q = i CA— z/A"^ < 6 < A , 

^ 2/i 6p+a^-2-a 1+^ 

CAT^6 if A 2(*P-«) < 6 < AR'H^, 

6fi+a/i-2-a v 

ARHi^A^PFirj < 6 

ifA^'^R^P = mm{A^/^R^^,R). 

Proof: From Proposition 13.21 we obtain that 

_ h^rx'^A-^^-f'^'^r^-hn) + g{A-^^-^^'^r^-^Xn) 

- y^(l-^)/2^1+/5 

is bounded from above. Moreover we know from the proof of Propo- 
sition |3]2] that h solves (with slightly reordered terms as compared to 
equation (|20ll ) 

b(rx' A-^i-^^/^j-i-iSy \ _ u(r) A-{i-fi)/2yi-py \ 

-2-^^ ^i-^^2/ + 2(fc(rx', A-(i-«/2ri-/^x„) • e„ + g')dnh 

\h{rx\A-^^-^^l^r^-Hn) + g'e„P |fc'(0, A-(i-^*)/Vi-^^x„)|2 



= -1-u-m-w -V - VI, 



where we have used the notation h' - (&i, &2/ •••/ ^n-i, 0) in the first line 
of the equation. 
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Let us first show that the right-hand is bounded: Using exactly the 
same argument as in the estimate of in Claim 3 of Proposition 13.21 
(consider the same argument with / := 1), we see that 

(40) sup \I\ < c{A'-3P+2p\m-2P) + Ai+/V2(«-/5)) < C. 
Next, 

\b{rx',A-('-Pyh'-Px„) - &(0,A-(^-^)/Vi-/^x„)| ,-„ 

(41) sup < sup j^^^ m\ 

< M^r"-^ sup|V/z(x',x,0|. 

Bi 

Similarly we may estimate 

(42) sup|ffl|< sup osc (-^^) < 4AV-2^ < C 

(43) sup \IV\ < 2A^^+fi^l^r"-^ osc Vj < C, 

(44) sup 1^1 < CA(^-«)+(^+")^'r^"-^^(i+«) < C, 

Bi 

(45) sup \V1\ < CA^^-m-»)y^(<^-fi-»fi) ^ Q 

Bi 

From go]), (|4l]), (gZl), dH]), (gH) anddHJ) as well as the parabolic 
comparison principle Lemma |6.4[ we obtain that |V/z| < C in Qi = 
{x; |x'| < 1, < x„ < 1}. 

Next we are going to estimate the oscillation of each term of the 
right-hand side in equation (|39|) in the x' variable: First we notice 
that the oscillation of VI in the x' directions is identically zero. 

Next, using the gradient bound of h, we obtain that for every 



< x„ < 1, 




(46) 


oscJ< CA^-^r^^ 


(47) 


oscJ/<2A(i+^')/V-^ 


(48) 


osc7JJ<4AV-^^ 


(49) 


osclV <CAr", 


(50) 


osc V < 2A^+^r2'«-^\ 
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So the oscillation in x' of the right-hand side in equation (j39|) is 
estimated hy c\ = C max , where we have used that in 

view of 2jS < « < 1, oscg^ I or oscg^ III is dominating when r is small. 

Now let w be the solution to 

IVwp -lWn = q + fixn) in Qi, 

w = h on ^Qi \ {x„ = 0}, 

where is chosen so that the right-hand side of the above equa- 
tion equals the right-hand side of (|39l) at \x'\ = 0. By the parabolic 
comparison principle Lemma [6T4l and the one-sided estimate Lemma 
16.21 for viscosity solutions we have for p in the super-differential of w 
at Xo that 

h{x) < w{x) < iu{xo) + p ■ {x' - x'q) + C\x' - Xgl^ 

< h{xo) + p-{x' - Xq) + C\x' - x'^f + Cq. 

Rescaling back to h'^ we see that this implies that (for some p in the 
super-differential of h'^) 

h\x',Xn) < h\Q,Xn) + p-x' + c(A^'-P^'^r'^l' q + A^'-P^'^r^-Vf) 

for < Xn < A~~r^~P and \x'\ < r. 

Using this estimate in an optimal way will yield the proposition. 
Rearranging terms and taking the supremum in B'^ on both sides 
leads to 

(51) snp(h^{x',x„)-h^{0,xn)-p-x') < dA^'-l'^'V^Pq+A^'-l'^'^rl'-'d^) 

provided that 6 G (0,r) and < x„ < A"t^6^"^. 

For fixed R and fixed 6 G (0, R) we want to find an r optimizing 

this estimate. Since we have a constraint that r < mm{A^^^R~^ ,R) it 
is natural to divide the proof into two cases. 

Case 1,R = mi{A^'^R^P ,R): That is, when 

R>A . 

When q = CA^r"~^^, or equivalently when 

(52) r < A"^, 
then we would want to choose 

P 2 

(53) r =A 2+a-2f!52+„-2/i_ 
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P 2 1-2/i 

(54) A 2+a-2p^2+a-2p <^ 4/i-« ^ 



From ((52]) we infer then that 
(54) 

or equivalently that 

(55) 6<A . 



Moreover we need r > 6 and r < R in order to use (pT]) , that is 

(56) 6 < A"^ 
and 

S 2+rt-2S 

(57) 5<AJR^. 

When (|52|) , (|56|) and (|57|) hold, then we can choose r according to (|53)) 
and deduce that 

sup(/z+(x',x„) - /z'^(0,x„) - p • x') < CA^"2(2+«-2ff5^'^2^. 

It still remains to consider the cases when at least one of (|52l) , (|56l) 
or (f57|) does not hold. 

If (|52l) holds but not (|57|) then we choose r = R and we deduce that 
— using q = Ah"~-^f^ which is equivalent to (f54)) — 

(58) sup (/z+(x', x„) - /z+(0, x„) - p • x') < c(A^R^+^'^^ + A^R^-^b^) 



^c(a-^R'^''-^+A^RP-'6^). 



In order to simplify (|58|) we use that 6 > A2R~^~ and thus 

A^R^-^62 > A^Ri+«-^ 
This way we may simplify (|58l) to 

(59) sup (/z+(x', x„) - /z+(0, x„)-p- x') < CA^Rl^-^6^. 



If (|52|) holds but (|56l) does not, then we choose r = 6 and deduce 
from dST]) that 

(60) sup {h^{x',Xn) - h^{0,x„) -p-x')< c(a^6^+«-^ + A^d^^^). 



The information that (|56l) does not hold implies 
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SO we may estimate 

We can thus simplify (|60ll to 

(61) sup (h*{x\x„) - h*{0,Xn) -p-x')< CA^6^+"-^. 

In case (|52l) does not hold — that is when q = CA^'f^r^^ — our first 
-JdL J_ 

choice of r is r = A 5 1+/5 ; notice that by the condition r < R, this is 
only possible for 

(62) 6 < A^Ri+^ 

Moreover we must have that r > 5 which in this case becomes 

(63) 6<A~W. 



So in case (|52|) does not hold but (|62|) and (|63|) do we obtain from (pTjl 
that 



(64) sup(/z+(x',x„)-h+(0,x,,)-p-^') ^ CA^^S^^ 

If neither ((52|) nor ((63l) hold then we chose r -2b and deduce that 

(65) sup (/z+(x',x„) - h^{Q,Xn) -p-x')< c(A^5^^^f^ + A^d^^f^). 

In order to simplify (|65l) we notice that as (|63ll does not hold, we have 
6 > A 2/* and thus 

A^6i-3^>A^6i-/^. 
Therefore we may write ((65ll as 

(66) sup(h+(x',x„)-/z+(0,x„)-p-x') < CA^6^+3/5. 

Finally, if neither ([52]) nor (|62]) hold, we use r = R in equation ((5T|) 
and 6 > A^i^^+^ to obtain 

sup(/i+-h+(0,x„)-p-x') < c(a^J^^+«-'' + A^J^''-162) < CA^6i+3^. 

This is the final estimate on Q in Case 1. 

Case 2, A^/^R^^ = mfiA^^^R"^ ,R): The argument is similar as in 
Case 1. In this case we have R < A~^. Again, if q = CA^r"~^^ then 
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the best balance in (jST) is obtained by the choice r = A 2+a-2pQ2+a-2is _ 
But when q = CA^r"'^^, then r < A~*F^' which is equivalent to 

6/i+a/i-2-rt 

(67) 6 < A 

when r = A 6 ^+"-^1' . In order to use ((5T|) we need r > 5 which with 
our choice of r is equivalent to 

(68) 6 < A"^. 

1 1 

Moreover, as r < A^R^-p , we must have 

2+,, 2+a~2p 

(69) 6 < A— R w 

for our choice of r. 

If (EZl, (Ell) and dUl) hold, then dHB implies that 



(70) sup(/z+ - /z"'(0,x„) -p-x')< CA^ 2{2«U/i)5i+2J-2/i_ 

As before we also need to investigate what happens if at least One 
of the conditions q = CAh'^~'^^, (|68|) or (|69|) does not hold. 

If (|67|) holds but (|69|l does not, then we choose r = A^R^-? and 
deduce from (jST) that 

(71) sup ih" - /z+(0, x„)-p- x') < c(A^ri+«-^ + A^R^^-^b^) 

= c(a— R— + R-^b^). 
On the other hand, the information that (|69l) does not hold implies 
that > A't^R—. We may therefore simplify ((71]) to 

(72) sup - /z+(0,x„) - p • x') < 



Next, if (|67| holds but (|68|) does not, we choose r = 26 and deduce 
that — using that (|68l) does not hold — 

(73) sup {h^ - h^{0,Xn) -p-x')< c(a^61+"-^ + A-^b^^l^) 



< CA— 

If (|67|) does not hold, then = CA^'h'^^ and the optimal choice of r 
in (|5T|) is given by 

(74) r = A 
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In order to satisfy 6 < r, we require furthermore that 

(75) 6<A"^. 

1 1 

And in order to satisfy r < A^R^-p we require that 

(76) 6 < AR^P. 

From (ITS) , (|76ll , (|74|) as well as the information that (|67l) does not hold 
we infer that 

(77) sup - h*{0,x„) -p-x')< CA^Pd^^^P. 

Notice that ([75]) is always satisfied in Case 2 as 6 < R < A ^p . 
We therefore only need to check what happens if neither (|76l) nor 
(|67|) hold in order to finish the proof. In this case we choose r = 

A5RT^. With this choice, (jSl]) implies that — using the fact that (|76l) 
does not hold in the last of the following inequalities — 

(78) sup (h^ - h^{0, Xn) -p-x')< C^A^r^+^z^ + j^l^^p-i^2\ 

□ 

5. Regularity for Obstacle Problems (Proof of Main Theorem) 

5.1. Heuristic arguments. In this section we combine the results 
from our previous sections and prove our main theorem, i.e. optimal 
regularity of the solution to the obstacle problem with Hamilton- 
Jacobi obstacle. 

The proof is again somewhat involved, so before we start let us 
describe the idea. 

Our goal is to prove optimal C^'^^^-regularity for minimizers to the 
double obstacle problem with h- as obstacles. The proof consists of 
several steps and an iteration argument. Here is a scheme of the 
steps: 

Step 1) From Corollary 12.51 and Lemma 16.51 we already have some 
regularity, in particular C^' 2^ -regularity for u and one-sided C^'^- 
estimates for h- at points where h- = u. 

Step 2) Having C^'^-estimates for u and one-sided C^'^-estimates for 
h-, we can apply Proposition |3]2] and Proposition 14. 11 which will give 
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US a slight, let us say an e, gain in the Holder exponent for one-sided 
estimates in the x'-directions for h-. That is, h- satisfy one-sided 
C^'^'^'^-estimates in the x'-directions. 

Step 3) Using that u solves the obstacle problem together with one- 
sided C^'^'^'^-estimates we can show that u G C^'^"^*^. In particular, we 
have gained an e = e{a,p) in regularity as compared to Corollary 12.51 
and Lemma [631 

Step 4) We can iterate Step 2) and 3) to gain more regularity of u, but 
in order to fully utilize the e gain in regularity we need to be able to 
control the C^'^^^-norm. of ii. As it turns out there is a constant E, such 
that if [fl]c« < ^ then we will get good control over the C^'^'^^-norm of 
u. So we will rescale u and a by a factor t, where t is chosen such 
that [fl(Tx)]cn < ^. With this rescaling we can iterate Step 2) and 3) to 
gain another e in the Holder exponent and also preserve the Holder 
norm. 

Iterating Step 2) and 3) we will see that u G C^'^ for any (i < a/2, 
but with a uniform bound on the C^'^-norm. The Theorem follows. 

In reality, the proof will be somewhat more involved as we have 
different regularity on different scales which will result in some tech- 
nical issues. 

With this strategy in mind let us turn to the proof of the Main 
Theorem. 

5.2. Proof of the Main Theorem. Without loss of generality we may 
assume that |Vm(0)| = |u(0)| = — Vu exists by Corollary 12. 5j — and 
that a{x) = e„ + b{x) with \b{x)\ < A\x\" in Bi. If this is not true, we may 
subtract u{0) + Vw(0) • x from u and /, add Vu(0) to a and rotate the 
coordinate system to obtain this situation. 

We may also rescale = ^^^withx = (^/A)^^" for a ^depending 
only on a and n to be determined later. The scaled solution 
will then minimize the Dirichlet energy in Br^ with Rj = 1/t and 
constraints /z* solving 

±\Vh^-a,\^ = ±1, 

where aj{x) = a{Tx). We set := [aj]o'(Bi/,) = ^■ 

From Corollary 12.51 we have oscg,. Uj < caI^^^'^"\^'^"^^^'^"'^ whenever 
r < A^^ mm{Rf*"'^^^, Rj) and also — using Lemma 16.51 — supg < 
CAy^^^^V^'^"/^^'^"^ That is enough to apply Lemma [6^ and Proposition 
SDwith |6o = ct/{2 + a) and r < mm{Al'^Rl'^^~^°\R,). 
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Let US assume, for now, that r < Rj < aJ . Then we have 

that for every x„ G {0,Aj ^ ri"/^o) and some p{x„) = p in the super- 
differential of h'^ at (0, x„), 

(79) sup Xn) - KiO, Xn)-p- x') < CaJ"'^'^^?'^^^ 



x'eB' 



if r < min(A, , A, , A/RT^ ). 

With 

F(r) = CA^'^'^r^^^Ts^^^ 

it follows from (|79l) that u satisfies the assumptions of Proposition 
12.41 that is, there is a constant C such that 

(80) osc,6B,.i/, < ca; 

if r < := min(Ai/2Ry('-^°U7^,A:^,A7Rr^). 

The estimate (|80|) gives us better regularity for small r. We would 
want to iterate that statement. Let us rewrite the statement with 

|6i = a/{2 + a- l^o): 

(81) oscb,m,< (CA, ^ )A,2 r^+^*i<A,^ r^+^i; 

here the last inequality is valid if A^ is small enough, however the size 
of At is independent of j6o and j6i (given that 1/2 > j6o, jSi > a/(2 + a)). 
As At = ^ by our rescaling, let us choose E, as the largest constant < 1 
such that the last inequality holds for all 1/2 > j6o, j6i > «/(2 + a). 
Then ^ depends only on a and n. 
We are going to apply Proposition 13.21 and then Proposition 14.11 

again with j6 = j6i. It is natural to split the proof into the four cases 

1 

(1) Rt > 1 (in which case A > ^) and < A,R';^\ 

(2) Rt < 1 (in which case A < 5) and R, > A^R^\ 

1 

(3) Rt > 1 (in which case A > 5) and R^ > A^Rl~^\ 

(4) Rt < 1 (in which case A < 5) and R, < A^R^° . 

Here (1) and (2) relate to the first and second block of values of Q in 

Proposition 14. 1[ respectively. These cases are stable in the following 

1 

sense: when (} increases in Case (1) or (2), then A^R^ ' increases if 
Rt > 1 and decreases if R^ < 1. That means that we stay in the same 
Case (1) or (2), respectively for larger j6. In Case (3) and (4) this is 
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no longer true, so we might have to use one of the blocks of Q for a 
finite number of iterations and then switch over to the other block. 
We will finish the proof in Case 1 first. 

Case (1) (R, > 1 and R, < Ay^R^^^"^"^): In this case we may iterate 
and get a (ii = «/(2 + a - 2jSi) such that (|80|) holds with (^2 replacing 
jSo- We may iterate indefinitely to obtain an increasing sequence of j6y 
such that fij = a/{2 + a - 2jSy_i), and (|80l) holds with fij replacing jSo 
for each ; G N. It is easy to see that (ij —> a /2 as /—> oo. Using that 
At = ^ is a constant depending only on a and n and that R^ > 1, it 
follows that 

(82) oscb,w, < CA^r^+'i 

for r < inf (min(A, , A, , AjR^^ )) < C„. 

Rescaling back to u we obtain from ((82]) as well as the definition of 
Tthat 

(83) OSCb,M = T OSCb,,,H, < CiT-"'V^"'^ < C2Al^V^'''\ 

_l 

for r < CsAj here Ci, C2 and C3 are constants depending only on a 
and n. 

_i 

When r > C^A^ " we obtain from jSo = that 

oscb,,m < C^A^^ r^+l^° < CsAy*^, 

where C4 and C5 are constants depending only on a and n. Combining 
the two estimates proves the Theorem in Case (1). 

Case (2) (R, < 1 and R, > a1^^rI'^^~^°\. 

WhenA^^Rf^ = inf(Ay^Rf^,KO/ we have to use the second block 
of values for Q in Proposition 14. II We deduce that 

sup \h*{x) - h*{0,x„) -p-x\< CA5"2(2+«-w/+2Ts% 
b;.(o,x„) 

ifr<min(A 2(¥o-«) , A A" R ^d-w J. 

Using Proposition 12.41 and iterating, we obtain that for f}j defined 
above, 

(84) osc.,B.H. < C^'"^/'^ 
for r < min(A7^, A:^, A^R;'"^'' ). 
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As in (|8T|I we conclude that 



(85) osC;,6b,"t <A J' ^ 



fSj 2+a- 
— 2+a 



for r < min(A7^,A,""'\A7^R;'""'' ). 

Sending ; to infinity and using that < 1 and that ^2{i-l^) 
decreasing in (ij we obtain that 

OSCx'eBU-i: ^ Q^T 

2+a 2 

for r < A^-* R^~" , where Ce is a constant depending only on a and n. 

Using Proposition l2.4l as well as the second case in the second block 
in Proposition 14. 1 1 in iteration, we also see that 

OSCx'eB'U-, < CyR'^r^ 

2+a 2 

for r > A^* Rt~" , where C/ is a constant depending only on a and n. 

Scaling back and using that At^ is a constant depending only on a 
and n we obtain that 



oscw < 



CsA^r^+i ifr<CioAf«, 
Cgr^ ifr>C«Af^, 



where Cg, Cg and Cio are constants depending only on a and n. This 
concludes the proof in Case (2). 

Case (3) {R, > 1 and R, > AI'^rI'^^~^°^): We start as in Case (2) and 
deduce that for all ; such that 



(86) ]^,>aJr"^-\ 
we have 



sup \h+{x) - h+{0,Xn) -p-x\<CA' 

6Pj+aPj-2-a Pj 2+a-2Pj 

forr<min(A '<*'^r«) , A"^, A^J^^^). 

If (|86l) holds for all ; then we are done as in Case (2). If (|86l) is 
not true then there is a largest /o such that the inequality holds for 

all < /o- It follows that R^ > \ and R^ < A\'^rI''''° . These are the 
assumptions in the iteration in Case (1), so we may proceed as in 
Case (1) and obtain the statement of the Theorem. 
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Case (4) (R, < 1 and < aI'^r]^^'^ We start as in Case (1) and 
deduce that as long as R^ < 

(87) osc...... < aJ"^/^-^, 

for r < mm(ArRl"'-^'\A, ^Ar'^AjRT^). 

If ^ A^^J^y*^ for all / then we are done. If there is a jo such 

that Rt ^ Ai^'^^R^^^ then we are in the situation of Case (2), so we 
may proceed as in Case (2) and obtain the statement of the Theorem. 
This finishes the proof. □ 



6. Appendix: Remarks on Viscosity Solutions For 
Hamilton-Jacobi Equations. 

In this appendix we will remind ourselves of some basic proper- 
ties of viscosity solutions for Hamilton-Jacobi equations. Most of the 
results in the appendix can be found in [?]. However the exposition 
in [?] is quite sketchy at points and some of the results that we need 
are not explicitly proved. We will therefore, for the readers conve- 
nience, supply some details, without claiming any originality. For 
the original and classical papers on the theory of viscosity solutions 
we direct the reader to [?], [?] and [?]. First we will state the definition 
of viscosity solutions for first order Hamilton-Jacobi equations. 

Definition 6.1. We say that u G C°(Q) is a viscosity solution to 

H{x, u, Vu) = f{x) in Q, 

if for every e Qandcp e C^(B,-(x°)) swc/i thatu{x^) = (p{x^) the following 
holds: 

(1) if u{x°) = (p{x°) and u{x) - (p{x) has a local maximum at x°, 
then H{x^,(p{x°),V(p{x°)) < f{x% 

(2) if t/(x°) = (^(x°) and u{x) - (p{x) has a local minimum at x^, 
then H(xO,(/)(x"), V(/)(xO)) > f{x°). 

This definition turns out to be the right one for first order Hamilton- 
Jacobi equations in the sense that it provides strong existence and 
uniqueness results. 

Viscosity solutions can be obtained by the vanishing viscosity 
method. That method also implies one-sided estimates for the second 
derivatives of solutions to smooth convex first order Hamilton-Jacobi 
equations which we will prove next. 



34 JOHN ANDERSSON, HENRIK SHAHGHOLIAN, AND GEORG WEISS 

Lemma 6.2. Let u e C°(Br x {-R, 0)) be a viscosity solution to 

^ + H{x, u, Vu) = m{x) 
ot 

such that H e and 

p ■ DpH{x,a,p) ■ p > c\pf- for every {x,a,p). 

Moreover we assume that m G C^'^. Then the distributional second deriva- 
tives satisfy 

^<^inB„j2X{-RM) 

i 

for all spatial directions Xj. The constant Q depends only on ||D^m||L«(Bj5)/ 
supg^ \Vu\ and H. 

A similar statement holds for the time independent case, that is if u e 
C^{Br) is a viscosity solution to 

H{x, u, Vm) = m{x) 

such that H e and 

p ■ DpH{x,a,p) ■ p > c\p\^for every {x,a,p) 
and m e C^'^, then 

d^u Co , 

for all spatial directions x,. The constant Q depends only on ||D^m||L"(Bj^), 
supg^ |Vm| and H. 

Proof: We will only prove the first statement for parabolic Hamilton- 
Jacobi equations. The proof for elliptic Hamilton-Jacobi equations is 
similar. Alternatively the elliptic proof may be derived by consider- 
ing f to be a dummy variable. 

Moreover, as the proof is well known and included only for the sake 
of completeness, we will consider the slightly simpler case H{x, a, p) = 
H{p). The general case can be handled similarly. 

Rescaling u{Rx,Rt)/R if necessary we may assume that R = 1. Let 
X e C°°{Bi X (-1,0)) be a non negative function such that x = ^ 
in Bi/2 X (-1/2,0) and x = close to {t = -1} and dB^ x (-1,1). 
Furthermore we may assume that |Vxl, |Ax|, |V;^p/|;t'l/l^/l^/<^^l < C in 
the support of x- 

We are going to argue by the method of vanishing viscosity. That 
is, we will approximate ^ + H{Vu) = m{x) by the equation 

du 

(88) - eAw + ^ + H(Vi/) = m{x) 

ot 



< — in Br/2 
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in order to deduce the desired estimate independent of e > 0. The 
Lemma follows from uniform convergence by letting e — > 0. 
Differentiating equation (|88l) twice in the direction we get 

-eAw, + ^ + H;{Vu) ■ Vuii + (Vui ■ H;^) ■ Vui = mu{x). 
Writing w = xuu we see that 

-xHp(Vi/) • Vuii + xmi - le^X ■ ^^ii - euuAx. 

At a point to) where w attains its positive supremum we have 
-Aw > 0, = and dw/dt > 0. The last two conditions are 
equivalent to 

VT / . X ""(^°, io)Vx(^°, ^o) 

= - 

Using this together with the convexity assumption on H, we end up 
with 

dx 

< Uii— — cx|Vm,P + UiiH'JVu) ■ Vx + XT^ii + 2e Uu - eiiuAx 

dt f X 

at (x°, io). Rearranging terms we get 

c;r|Vu,f < + H;(ViO • Vx + 2e^ - eAx)M, + m,x. 

Using that — ^by our choice of x arid assumptions on H — the terms 
in the parenthesis may be estimated by a constant C depending only 
on supg^ |Vw| and H, and observing that the final term mux is also 
bounded by the assumption m G C^'^, we obtain that at the point 
where w attains its supremum, 

cxNuif < Quiil + C, 

where C depends on supg^^^.^^ m,,, supg^^^.^^ |Vw| and sup |Hp(Vt/)|. 
Multiplying both sides by x iiriplies that \iv\ < C. Which in terms of 
Uii becomes 

sup Uii < sup xuii < C. 

Bi/2X(-1A0) Bix(-1,0) 

□ 

Later on we will also need good stability estimates for solutions, 
proved in the next two Lemmas. 
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Lemma 6.3. Let H e and let u G C°'^(Br) and v G C'''^(Br) he viscosity 
solutions in Br to 

H{x, Vu) = n{x) 

and 

H{x,Vv) = m{x) 

where 

f ■ D^H{x,p) ■ p > c\pffor every {x,p), 
m > n > A > 0, u < V on dB^ and ||m||l"(Bj,), II^^IIl~(Br) ^ K. Then 
u <v <u + C(A, K)R max ( sup(m - n), sup{v - u)). 

Br dBR 

Proof: It is sufficient to prove the Lemma when sup^g^(m - n) and 
supg^(u - u) are small. By a rescaling we may assume that R - 1. 
Claim 1: Let 

= e^+i - e-\ 

Then there exists a function F{x, a, p), convex in p, such that 

dF(x, a, p) 

(89) — - ^' 

and in the viscosity sense 

F{x, Wi, Vwi) = 0. 
Moreover Wi = e^*^ - e~^ is a viscosity solution to 

F{x, 102,^102} = (e^'^^ - zy2)(w7 - n). 
Proof of Claim 1: A simple calculation shows that 

(e^+i - w^)H{x,Vw^/{e^+^ - w^)) = {e^^^ - Wi)m, 
so if we denote 

F{x,a,p) = {e^^' - a)H{x, - (e''^' - a)m, 

the first statement in the claim follows. That F is convex in p follows 
from the convexity of H. A direct calculation yields 

where we have used convexity of H. That 

f(x,W2/ Vwi) = {^^^^ - W7)j{m - n). 
follows from a simple calculation. The claim follows. 
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Next we observe that we may assume that Azfi < Ca for some 
constant Ca: we may regularize f so it becomes C^. Denote the 
regularized version of F by f ^ and let wl solve 

F^{x,iol,Vwl) = 0, 

then by Lemma [6^ Aw^ < C^. We also know that e C" uniformly 
for all T > 0, thus we find a subsequence Tj — > such that w^^' — > zf ° 
uniformly Using uniqueness results for viscosity solutions it is easy 
to see that = Wi. So if we can show the Lemma for each zf^' the 
result follows for Wi. Thus we may assume that Awi < Ca as long 
as our final estimate does not depend on Ca- We will make several 
more regularizations in what follows. In order to simplify notation, 
we will by the just explained argument assume that Azf i < Ca, and 
continue to work with zfi and not with zf^'. 

Next, we apply a standard convolution type regularization to Wi 
and u>2, that is we denote 

= iOi{y)(p6{x - y)dy, 

where is a standard moUifier. Then 

F{x,w\,Vw\) = hi 

and f (x, zy^, Vzy^) = (e^+^ - w^{m -n) + h^, 

where — > uniformly in U for all p 6 [1, +oo) as 6 — > 0. To do this 
regularization we need to assume that u and v are solutions in B1+5 
or prove the result in Bi-s. However as 6 — > this will not make any 
difference so we may as well ignore this slight complication. 
Finally we will denote by w^^ the solution to 

-eAzyf + f (x, zyf , Vzyf ) = h^2, zyf = zy^ on dB^. 

With these regularizations we end up with 

-eA(z(;f - wl) + f (x, zyf , Vzyf ) - f (x, zy^, Vzy^) 

< {e^*^ - w^{m -n) + cCa + hl- h^. 
Next, in order to linearize we define 

B{x)= f F,j,{x,wl,tVw^^' -{l-t)Vzvl)dt. 
Jo 

Then, taking equation (|89l) into consideration, 

-eA(z(;f - wl) + B{x) ■ (Vzyf - Vzy^) + A(z(;f - wl) 
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< - iU2){m -n) + cCa + h^^- h^. 

Now let /^'^ be the solution to 

-eAf + B{x) ■ + Af'' = - W2){ni -n) + CaC 

with boundary values p'*^ = sup^g {lu'^'^ - w^'*^) and g^'^ be the solution 
to 

-eAg^'' + B{x) ■ Vg^'' + Ag^'' = h\- h\, 
with boundary values g^'^ = on dBi. By the comparison principle 

(90) - < f'' + g^''. 

Since - — > in LP for any p < +oo as 6 — > we see that 
supg^ \g^'^\ — > as 6 — > 0. Next we notice that where attains its 
supremum we have \Vf^'^\ = and Af'^'^ < 0. It follows that 

sup/^''^ < max ( sup /^'^, sup yfe^"^^ - iV2){ni - n) + —). 

Letting first 6 — > and then e — > we see that equation (|90l) implies 
that 

W2-ZV1 < max ( sup(zy2 - Wi), sup — fe^"^^ - 'W2){ni - n)). 

b, ' ' 

Writing this last inequality in terms of u and v we get 

e"" - e"'" < max i sup (e"" - e~^\ sup -^(^ - "))• 

Thus 

1 - e"-" < jj^ax ( sup (1 - e"""), sup -(m - n)). 

^ 5Bi ^ Bi A ^ 

As pointed out in the beginning of the proof it is enough to show 
the Lemma when supg^(m - n) and sup^g^(z; - u) are small. Using 
< 1 - < 2E, for small E,, we end up with 

sup(u - u) < 2e^^ max f — sup(m - n), 2 sup(?7 - u)j. 

□ 

Next we need a comparison estimate for parabolic Hamilton-Jacobi 
equations. The proof is very similar to the proof of Lemma [63] so we 
will omit it here with a reference to [?]. 
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Lemma 6.4. Let H e and let u e C^-^Br x {-R, 0)) and v e C°'HBr x 
{-R, 0)) be viscosity solutions in Br x {-R, 0) to 

+ H{Vu) = m{x) 

ot 

and 

^ + H{Vv) = n{x), 

where 

p ■ D^H{p) ■ p > c\p\^for every p, 

n > m, u = V on dBRfor t e {-R,0) and u{x, -R) = v{x, -R) for x e Br. 
Then 

u < V < u + CR sup (n - m). 

BrX(-R2,0) 

Proof: For a discussion of a proof see [?]. □ 
From the regularity results in Lemma 16.21 and 16.31 we can easily 
deduce some elementary, non-optimal, one-sided estimates for so- 
lutions to Hamilton-Jacobi equations, even when the data is not C^. 
This incidentally provides our starting regularity for a bootstrap ar- 
gument which in turn yields optimal regularity. 

Lemma 6.5. Let h be a viscosity solution to \Vh - a\^ = 1 in Br, a G 

C"(Br;R") and [(i{x)]Qa(g^') < A. Then for any x° G Br/2 the super- 
differential ofh at x^ is not empty, and for any p in the super-differential of 
hat x^ we have 

(1) 

sup [h{x) - p • (x - x°) - h{x^)) < CA^I(^+")/+»li^+») 

for r < min{A^/^R^^+''^/^,R), 
(2) and 

sup (h{x) - p • (x - x°) - h{x^)) < d— 

forA-^l" > r > Ai/2j^(2+a)/2_ 
The constant C depends only on n. 
Proof: We may assume that h{0) = 0. 

First we notice that —as r < R— if r > A'^/'' then r < A^/^R^^+"^'^. 
Also, if r > A-^'" then 

1 = sup \Vh - a\> sup \Vh\ - sup \a\ = sup |Vh| - Ar". 

Bj- Bf Bf By 
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Using that Ar" > 1, it follows that 

sup \h\ < r sup \Vh\ < r(l + aA < CA^/i2+"^r^+'^/(2+«)_ 

By By 

This shows that if r > A~^^" then both the assumption and the con- 
clusion in (1) hold. 

We therefore only have to show the Lemma when r < A~^^". We 
do this by a barrier type argument. We may assume that a{0) = 0. 
We also define the barrier w as the solution to 

|Vw|2 = 1 + {2A6" + A^6^") in B^ix"), 
w = h on dB^ix^), 

where 6 is to be determined later. Then \Vw\ < ^^/T+TmS^^Ta^S^") < 
1 + (lAd" + A^6^") and thus 

\Vw -a\> \Ww\ - \a\ > Vl^r(2A6«TA262^ - Ad" 
= 1 + ( Vl + (2A6« + A262«) - (1 + A6")) = 1, 
so w is a super-solution to \Vh - = 1. Similarly, 
\Viu -a\^<l + CA6"[l + A6"). 

We may thus use the comparison estimate for Hamilton- Jacobi equa- 
tions (Lemma 16. 3|) and derive 

h<w< w(x°) + p ■{x-x^)+ ^ 



< h{x°) + p • (x - x°) + ~ + CA6^^"(l + Ad") 

for p in the super-differential of w. Rearranging the terms and taking 
the supremum in Br{x^) for some r < 6 yields 

2 

(91) sup (h{x) - p • (x - x°) - h{x^)) < cC- + A6^^"(l + A6")). 

We want to find the right balance between r and 6 that optimises ((9T|) . 
It is convenient to divide the last part of the proof into two cases: 

Case 1: When r < min(A^^^R^, A"n), we use 6 = A~2^'r~^ in 
(HD and deduce that — using r < A"^^" which implies that A6" = 

/\_2+af2+a < \ 

sup {h{x) - p • (x - x°) - h{x^)^ 

B,.(xO) 

< c(A5ijri+2fs + A6^+"(l + A6")) < CA^^r'^*^^. 
From the definition of r and 6 it is easy to check that r < 5 < R. 
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Case 2: Next, if A^'^R-^ <r< A'^l" then we use 6 = R in (HD and 
deduce that— using A ^/^j^T^ < r,A^I^R^ < r < A"!/" and AR« < 1— 



2 2 

sup {h{x) - p • (x - x°) - < C(^ + ARi+«(l + AR«)) < C^. 

Observe that this interval is empty unless R < A~^l". □ 
We end this appendix with a lemma reminiscent of a Whitney 
extension theorem, which will be useful in the main text. 

Lemma 6.6. Let h G C°'^(Bi) and assume that h satisfies the one-sided 
C^'^ -estimate 

(92) h{x^ + x)< h{x°) + p,o ■ {x - x°) + Qlxr^ 

for every p^o in the super-differential ofh at x° and every x° G Bi. Moreover 
assume that h{x) > -Ci|xp'^'^ and that h{0) = 0. Then 

(93) |p,| < E(j6)(Co + Ci)|x|^ /or x G By2; 
here E(j6) depends continuously on pfor p G (0, 1). 

Proof: Take y G B,- and let py = p he in the super-differential of ?z at 
y. As is semi-concave by (|92l), we know that the super-differential 
of is non-empty for every x°. Notice that, due to (|93l) , the super- 
differential of at the origin contains only the zero vector. 

Therefore h{y) < CqT^'^^. Now consider z = y - ep with 

e = (6|pp-f 

for some small constant 6. By one-sided estimates from above and 
below we have 

< -Cily - ep\'^^ < h{z) < h{y) - e|pp + Qe^+^IpP^^ 

where is a constant depending only on j8 and n. Reordering terms 
and using that h{y) < Cor^+^ yields 

[b'l^ - (CpC, + Co)6^'^p'^^)\p\^ < (Co + q,c^y^^. 

Choosing 6 = (2Co -i- 2C^Ci)"^ the previous inequality becomes, for 
some C depending only on j6, 

< c(Ci + Cof^l'^'^r'^^. 
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Taking both sides to the power of j6/(l + (}), the lemma follows. □ 
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